A PHASE-FIELD APPROXIMATION OF THE WILLMORE FLOW 
WITH VOLUME AND AREA CONSTRAINTS 



PIERLUIGI COLLI AND PHILIPPE LAURENgOT 

Abstract. The well-posedness of a phase-field approximation to the Willmore flow with area and 
volume constraints is established when the functional approximating the area has no critical point 
satisfying the two constraints. The existence proof relies on the underlying gradient flow structure of 
the problem: the time discrete approximation is solved by a variational minimization principle. The 
main difficulty stems from the nonlinearity of the area constraint. 



1. Introduction 

Biological cell membranes define the border between the interior of the cell and its surrounding 
medium and can be roughly described as a lipid bilayer in which several kinds of lipids are assembled 
and through which proteins diffuse. The size of the cell (a few microns) is typically much larger than 
the thickness of the membrane (a few nanometers) and a possible approach to model the geometric 
properties of the latter is to assume the membrane to be a two-dimensional embedded surface S in 
with a shape at equilibrium being determined by the Canham-Helfrich elastic bending energy 
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see, e.g., [H O HH [13] and the references therein. Here, % := {hi + h2)/2 is the arithmetic mean 
of the principal curvatures hi and h2 of S (scalar mean curvature), /C := hih2 is the product of its 
principal curvatures (Gaufi curvature), k and kg are the bending rigidity and the Gaussian curvature 
rigidity, respectively, and T-Lq denotes the spontaneous curvature which accounts for the asymmetry 
of the membrane. Let us mention here that, when I-Lq = kg = and k = 1, the functional Sch is 
nothing but the Willmore functional which is a well-known object in differential geometry [5T]. Two 
natural geometric constraints come along with cell membranes: the inextensibility of the membrane 
fixes the total area while a volume constraint follows from its permeability properties O [T2] . 

Recently, experimental results have shown evidence of dynamic instabilities in membranes, see, e.g., 
[1] and the references therein, and provided the impetus for the development of dynamical models. 
A first approach is to consider the gradient flow associated to the Canham-Helfrich functional which 
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describes the time evolution of a family of (smooth) surfaces {T,{t))t>o and reads (in the simplified 
situation k = 1 and T-Lq = kg = 0) 

V = -AsH -2n{n^ -K) , (1.2) 

where V and denote the normal velocity to S and the Laplace-Beltrami operator on S, respec- 
tively. Volume and area constraints can also be included and result in an additional term in the 
right-hand side of (11.21) of the form i + rnH, the parameters I and m being the Lagrange multipliers 
associated with the two constraints. The main drawback of this approach is that it requires to solve 
a highly nonlinear free boundary problem which is difficult to study analytically and costly to com- 
pute numerically. However, numerical schemes have been recently developed for geometric evolution 
equations such as (II. 2p . see, e.g., [2]. 

A well-known alternative to free boundary problems is the phase-field approach where the sharp 
interface E(t) is replaced by a diffuse interface which is nothing but a thin neighbourhood of thickness 
e of the zero level set of an ^-dependent smooth function, the order parameter. For biological 
membranes, this approach has been developed in several recent papers with and without the volume 
and area constraints [3l HI [5], [6l [11], [121 [E] and can be described as follows, still in the simplified 
situation k = 1 and T-Lq = kg = 0: let be a smooth double-well potential (for instance, W{r) = 
(1 — r^)^/4) and for e > and w G H'^{Q) define the free energy by 

Es[w]:='- dx, (1.3) 

where Q is the spatial domain which comprises the cells and their surrounding medium. The corre- 
sponding phase-field model is the gradient fiow of in and reads 

d,v, = £ A/i, - - W"{v,) /i, , fi, = -Av, + \ W'{v,) , (1.4) 

supplemented with suitable initial and boundary conditions. As e approaches 0, the function is 
close to the values ±1 in large regions of the domain separated by narrow transition layers of width 
e around the zero level set {x : v^it^ x) = 0} of at time t. It is this time-dependent family of level 
sets which is expected to converge as e — to a family of surfaces (T,(t))t>o evolving according to the 
geometric motion (11. 2p . Formal asymptotic expansions have been performed to check the consistency 
of the free energy (II. 3p and the phase-field model (11.40 with (II. ip and (II. 2p . respectively, in the limit 
e = [5l [TOl [HI [20]. While no rigorous justification seems to be available so far for the evolution 
problem, the relationship between the minimizers of the Canham-Helfrich functional (ll.ip (without 
constraints) and those of the free energy (ll.Sp has been the subject of recent studies [3l [9l [161 [18] . 
In fact, the existence of minimizers of the free energy (II. 3p follows by standard arguments from the 
theory of the calculus of variations. But the well-posedness of the evolution phase-field model (ll.4p 
for a fixed positive e seems to be less obvious and, as far as we know, it is a widely open topic, 
though the phase-field approach has been used quite extensively to perform numerical simulations 
of the dynamics of biological membranes [51 [HI [HI [12] • We are only aware of two contributions in 
that direction. On the one hand, the well-posedness of the phase-field model (II. 4p with e = 1 and a 
volume constraint fixing the average of v has been shown in [8]. On the other hand, the existence of a 
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weak solution to a system coupling a phase-field model (similar to ( ]1.4p but with a convection term) 
with the Navier-Stokes equation is established in [22], relaxing the volume and area constraints by a 
penalisation approach. Therefore, accounting for both volume and area constraints in the phase-field 
approach to biological membranes does not seem to have been considered yet and is the focus of this 
paper. Before describing precisely our result, we recall that, in the phase-field approximation, the 
volume and area conservations read O [10] : 



Vsit,x) dx = const, and -Fe[t'e] •= / ( ^ \^'"e{t,x)\'^ H — W{vs{t,x)) ] dx = const. (1.5) 



Indeed, recall that, as e — )■ 0, the functional approximates the perimeter functional 

The purpose of this paper is then to investigate the existence and uniqueness of the phase-field 
approximation to the geometric flow (11. 2p when both volume and area are fixed as described in (II. 5p . 
We may then set e = 1 in the forthcoming analysis and define 

W{r) ■=^{r^ reM, 

F[w]:=- W{w{x))dx, weH\n), (1.6) 

2 Jn 

where a is a given positive real number. The phase-field approximation of (II. 2p with fixed volume 
and area turns out to be 

dtv = Afi-W"{v) fx + A + B fi, (t, x) G (0, oo) X 1] , (1.7) 

fx = -Av + W'{v), (t,x) e (0,oo) X (1.8) 

■ n = V/i ■ n = 0, {t,x) e {0,oo) X dft, (1.9) 

v{0)=vo, xeQ, (1.10) 

where n denotes the outward unit normal vector field to dfl and W the spatial average oi w E L^{^), 
that is, 

w := — — / w{x) dx . 



In (11. 7p . A and B are time- dependent functions and the Lagrange multipliers corresponding to the 
volume and area constraints 



v{t)=v^ and F[v{t)\ = F[vq] , t>0. (1.11) 

When the area constraint is not taken into account (which corresponds to take i? = in (II. 7p and 
keep only the first constraint in (II. lip ), the well-posedness of the resulting version of (ll.7p - (ll.lip is 
shown in [8]. In that case, the equation turns out to be a gradient flow for the free energy 

E[w]:=- [ [-Aw + W'{w)f dx, weH\n), (1.12) 
2 Jn 

and the existence proof exploits this structure and relies on a time-discrete minimization scheme. 
This gradient flow structure is still available for (ll.7p -f im|) with the main difference that we now have 
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two constraints including the additional one, which is nonlinear and generates several new difficulties 
in the analysis of the minimizing scheme. In particular, we emphasize that the two constraints may 
not be linearly independent and this happens in particular for critical points of F under a volume 
constraint. Deriving the Euler-Lagrange equation for the time-discrete minimization scheme is then 
not obvious, with the further drawback that the area constraint is nonlinear. This difficulty strikes 
back when we wish to estimate the Lagrange multiplier B and we have to restrict our analysis to the 
case where critical points of F under a volume constraint cannot be reached during time evolution. 

Let us now introduce some notations: given a G M and (5 G [0,oo), the fact that there exists at 
least one function w G H^{VL) satisfying simultaneously w = a and F[w\ = /3 is not granted and 
requires a compatibility condition which we describe now. We set 

/3„:=inf{FH : w E H\n) , w = a] , (1.13) 

which is well-defined owing to the nonnegativity of F, and 

Ml^p := [w G H^in) : w = a and F[w] = /?} , (1.14) 

:= {w e Hf^in) : weMl^} , (1.15) 

where 

Hlj{n) := {w G H\n) : Vv-n = on dQ} . 
Clearly, A^^^^ = if /3 G [0,/3q), while we shall show below that -Mq^/j, ^ = 1,2, is quite large if 
(3 > f3a (see Lemma [2?T] below) . 

Finally, as already mentioned, we have to exclude some values of the parameters {a, /3) for which 
there are critical points of F under a volume constraint in A^^^. To this end, we introduce the set 
Za^p defined by 

Z^^^ := ^^w e Ml^/s ■ -Aw + W'{w) - W'{w) = in fi} , (1.16) 

and shall require this set to be empty, an assumption which is fulffiled if /3 is sufficiently large 
compared to see Lemma [2.31 below. We may now state our result: 

Theorem 1.1. Consider a G M and /3 G (/3q,, oo) such that 

^a,/3 = 0. (1.17) 

Given an initial condition vq G A^^^, there is a unique function 

V G C([0, oo) X H) n L°°(0, oo; H\n)) , v{0) = Vq , 

such that, for all t > 0, 

v{t) G , /i := -Av + W'iv) G L\0, t; Hl,{n)) , 

and there are two functions A G L^(0,t) and B G L^(0,t) such that 

dtv = Aii-W'\v) 11 + A + B II a.e. m (0,t)xl]. (1.18) 

In fact, there hold 

A + BJI= W"{v)fi (1.19) 
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and 



B — ji 



I = + [ W"{v) dx - W"{v)fx [ fxdx 



(1.20) 





s e [0, t] . 



(1.21) 



Observe that the identity f ll.20p defining B is meaningless if /i(to) is a constant at some time to, 
that is, if w(to) ^ ^a,/3- The main purpose of the assumption (11.171) is then to prevent this situation 
to occur. 

A further consequence of our analysis is the time monotonicity of the free energy along the fiow 
which is a natural outcome of the gradient fiow structure of (ll.7p - (ll.lip . 

Corollary 1.2. Under the assumptions and notations of Theorem \l.l[ the map 



The outline of the paper is as follows: in the next section, we collect some preliminary results 
concerning the structure of A^^^, i = 1,2, some functional lower bounds on F and E, and the 
fact that Za,i3 is indeed empty for f3 large enough, along with a useful functional inequality in that 
case. Section [3] describes the minimizing scheme for one time step. Estimates are also derived there, 
allowing us to pass to the limit as the time step decreases to zero and obtain the existence part of 
Theorem 11.11 in Section |H The uniqueness part of Theorem 11.11 is proved in Section |5] and the proof 
heavily relies on the positivity property (ll.2ip which allows us to control the difference between the 
area Lagrange multipliers of the two solutions. 

2. Preliminaries 

Let us first show that -Ma^p, "^ = 1,2, is quite large when a G M and (3 G {/3a, oo) as claimed in the 
Introduction. 

Lemma 2.1. Consider a G M. There is at least a function Wa G H^f^Q) such that = a and 
F[wa] = Pa- In addition, if [3 E (/3q,,oo) and ip G H^iVt) (resp. if G HlfiVL)) satisfies ^ = and 
(y9 ^ 0, then there zs A > depending on ip such that Wq, + A<y5 G M.]^^^ (resp. Wq, + A</9 G M-I^^p)- 

Proof. Let a G M. The existence of Wa follows from the nonnegativity and weak lower semicontinuity 
of F by classical arguments of the theory of the calculus of variations; moreover, Wa solves the 
Euler-Lagrange variational identity associated with F, that is. 



whence Wa G H1j{VL). Consider next [3 G {13a, oo) and a function ip in either H^{n) or HI{Vl) such 
that = and ^ 0. Introducing the function / defined by /(A) := F[wa + Ay?] for A > 0, we 
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realize that / is a continuous function in [0, oo) with /(O) = = /3q, < /3 and, thanks to the 

nonnegativity of W , 

/(A) > \ ( ||Vw„||^ + A^llVv^lls + 2A / -V^ dx] — ^ oo . 
1 \ Jn J 

The mean-value theorem then guarantees that there is at least A^ > such that /(A^) = (3, that is, 
Wa + A<^(y9 G -M-a^is for either i = 1 or i = 2. □ 

We next show that the functionals F and E control the if^-norm and the if^-norm, respectively. 

Lemma 2.2. Given a G M, there is Ci > such that 

\m <Ci(^l + v^FH) for w G H\Q) such that w = a, (2.1) 

\h2 <Ci(^l + ^/E[w]^ for w G Hfj{n) such that w = a. (2.2) 

Proof. Let w G H^{Q) with W = a. We readily infer from the Poincare-Wirtinger inequality 

llty — W||2 < C2 II Vw||2 , (2.3) 

and the nonnegativity of W that 

Ikll^i < 2 {\\w-W\\l + ||in||2) + ||Vu;||^ < (1 + 2C2) \\Vw\\l + 2a^\Q\ 

< (2 + 4C2) F[w]+2a^\n\ , 

whence (12.11) follows. Next, let w G Hff{Q) with w = a and set /i := —Aw + W'{w). The definition 
(|LT2|) of E entails ||/i||^ = 2E[w] so that /i G L'^ifl) and 

||u; — a||2 ||/i||2 > {w — a) dx = || Vty||2 + a {w — a) {w^ — w) dx > || Vtf II2 — C . 
Jn Jn 

Thanks to (12.31) and Young's inequality, we further obtain 

II Vw||^ <C + C2 \\Vw\\2 Wi^h < + C (1 + , 

and consequently 

||Vw||^ < C{l + E[w]) . 

Using once more (12. 3p gives 

\\w\\l, < C (1 + E[w]) . (2.4) 
Finally, observing that r i— )■ W'{r) + ar is non- decreasing, it follows from the definition of fi that w 
solves 

—Aw + W'{w) + aw = fi + aw in Q 

with homogeneous Neumann boundary conditions, and a classical monotonicity argument ensures 
that 

II Aw II 2 < WfJ- + aw\\2 < \\fi\\2 + a \\w\\2. 
Combining this estimate with (12. 4p readily gives (12.21) . □ 
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The last result of this section is devoted to the set Z^^p, defined in (I1.16p . We prove in particular 
another fact claimed in the Introduction, namely that, given a G M, the set Z^^ is empty at least 
for /3 large enough, so that Theorem 11.11 can be applied in that case. Throughout the paper, we use 
the following notation: given w G satisfying W = 0, the function Af{w) G H'^{Q) is the unique 

solution to 



- AJV{w) =w in Q, VAf{w) ■ n = on dQ , satisfying Af{w) = . (2.5) 

Lemma 2.3. Consider a G M and set Oa '■= {P G {f3a, oo) : Za,p = 0}- 

(1) The set is open and /3 G Oa if P is large enough. 

(2) Assume that /3 G O^- Given M > 0, it turns out that 

2 



niM '■= inf 



VA/" i^-Av + W'{v) - W'{v] 



: veMis, E[v]<M}>0. 



(2.6) 



Proof. (1) Let us first show that (/3q, oo) \ Oa is closed. Consider a sequence (/3n)n>i in (/?«, oo) such 
that /3„ ^ Oa for each n > 1 and /3„ — )■ /3 G (/?„, oo) as — )■ oo. Then, according to the definition 
of Oa, for each n > 1, there is Wn G A^a,/3„ such that —Awn + W'{wn) — W'{wn) = in il. Since 
{Pn)n>i is bounded, it follows from Lemma [2^ that (wn)n>i is bounded in if^(r2). Furthermore, the 
properties of Wn and the continuous embedding of H^{Q) in L^{Q) entail that 



EK] = 2 iy'(?i;„) < C (1 + \\Wn\\t) <C{1+ \\Wn\\%i) , 



so that {E[wn])n>i is bounded. A further application of Lemma 12.21 ensures that (wn)n>i is bounded 
in H'^[Q). We then deduce from the compactness of the embedding of Hjf{Q) in H^{Q) fl C{Q) that 
there are w G Hjf{Q) and a subsequence (ifnfe)fc>i of (w„)„>i such that 

in H^{n) and C(n) , ^ w in H^{n) . 



It is then straightforward to check that w G M-a fs satisfies —Aw + W'{w) — W'{w) 
that is, w G Za^/s- Thus, /3 ^ Cq. 

Consider next v G -2a,/3- Then, on the one hand, we have 

/ V {-Av + W'{v)) dx = \\Vv\\l + a v^v'^ - l) dx = 2F{v) + - {v^ - 1) dx 

= 2P + ^{\\v\\t-\n\) . 

On the other hand, it results that 

/ V {-Av + W'{v)) dx = W'{v) V dx = aa {v^ - v) dx . 
Jn Jn Jn 

The above two identities give 

2i3 + ^{2a'-l)+- [ {v^-2av') dx = 0. 
2 Jn 



in Q, 
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Since — 2ar^ > — 27a'^/16 for r G M, we deduce that 

^ a\n\ / „ 27 A 

2 16 ) -°' 

which is not possible if P is large enough. 

(2) Assume for contradiction that rriM = and let {vk)k>i be a minimizing sequence. Since 
Vk e -^0,13 with E[vk] < M, it follows from Lemma [2l2] that {vk)k>i is bounded in Hf^^Q) and thus 
compact in H^{Q) and C{Q). Therefore, there are v G Hjf{Q) and a subsequence of {vk)k>i (not 
relabeled) such that 

Vk — >v in H\Q) and , ^ w in H^{Q) . 

These convergences readily imply that v = a and F[v] = /3, so that v G ■'■'^ addition, setting 

ffc := — At>fc + Vr'(wfc) — W{vk) for > 1, we also deduce that {i'k)k>i converges weakly in 
towards := —Av + W'{v) — W'{v) while the property niM = entails that {i'k)k>i converges to zero 
in H^{Q)'. Consequently, 1^ = 0, from which we conclude that v G Za,i3 and get a contradiction. □ 

3. The minimizing scheme 

We fix a G M and /3 G (/3„, oo). Given / G L'^{Q) and r G (0, 1), we introduce the functional Jvj 
on Hj^lQ) defined by 

\\w - fP 

J'rAw]:=^-^^^ + T E[w], weHfjin), (3.1) 

and consider the following minimization problem 

UJrj:=mi{j^rj[w] : w e Ml^} . (3.2) 

Since E is non-negative and Ai"^ ^ is non-empty by Lemma [2m Urj is well-defined and non-negative. 

Lemma 3.1. The functional Trj has at least a minimizer in M^p- ^''^ addition, any minimizer v 
of Ft J in M-I^p satisfies 

E[v] < E[f] . (3.3) 
Proof. For each k > 1, there is G A^^^ such that 

COrJ < J'rAVk] = + r E[Vk] < UJrJ + \ ■ (3.4) 

We introduce fik := —Avk + W'{vk) for k > 1. From Lemma [272] and (13. 4p it follows that {vk)k>i is 
bounded in Hjf{Q). Since H^{Q) is compactly embedded in H^{Q) and C{fl), there are v G i/^(f2) 
and a subsequence of {vk)k>i (not relabeled) such that 

Vk — y V in H^(Q) and C(r2) , ^ f in iJ^(f2) . 



A PHASE-FIELD APPROXIMATION OF THE WiLLMORE FLOW WITH VOLUME AND AREA CONSTRAINTS 



9 



These convergences readily imply that v = a, F[v\ = (3, and /i^ ^ /i := —Av + W'{v) in L'^{Q). 

42 



Consequently, v G A1a,,3 

Wr,/ < -^rjH = r,'' '"^ + 2-r ||/i||2 < lim + liminf 2r ||/ifc|i2 



\^-f\\l , o_n..n2/ IK -/IP 



2^ 

fc— >oo 1^ fc— >oo 

liminf fe] = cu^,/ - 



We have thus established that w is a minimizer of Jvj in □ 

Remark 3.2. A similar argument shows that the free energy E has at least a minimizer in A^^^. 
Constructing a minimizer to the Canham-H elf rich functional f ll.ip turns out to he far more com- 
plicated even without constraints (see [T71 [19] and references therein). Existence of axisymmetric 
minimizers to this functional with volume and area constraints has been recently proved in . 

The next step is to derive the Euler-Lagrange equation corresponding to the minimization problem 
(13. 2p . At this point, a new difficulty shows up as the constraints are not always independent. Indeed, 
setting I[w\ := W, the differentials of / and E are DI[w\ = and DE[w] = —Aw + W'{w) and 

are not linearly independent if —Aw + W'{w) is a constant. Assuming that this is not the case, we 
have the following result: 

Lemma 3.3. Assume that v G M\ p solves the minimization problem fl3.2p and is such that fi : = 
—Av + W'{v) is not a constant. Then fi G if^(fi) and there are real numbers A and B such that 

V - f 

- Afi + W"{v) fi = A + B fi m Q. (3.5) 

T 

Proof. Owing to the nonlinear constraint E[w] = (3, we proceed as in [231 Proposition 43.6]. We 
observe that, since fi is not a constant, 

-71)2 = ll/i- /ills > 0. (3.6) 

Next, it turns out that jJ^-Jt is not sufficiently smooth for the forthcoming analysis and a regularization 
is needed. Owing to the density of C^{fl) in L'^{fl), for 77 G (0, 1), there exists i/^ G C^{^) such that 

1^ = and ||z/^ - (/i - 71)112 < ?7 • (3.7) 

Combining (13.60 and (13.70 gives 

(/w, i^v)'2 = (l^-T^^ ^v)2 = ll/u - 71II2 + (/i - 71, ly^- {fi- JI))2 > ll/i - 71II2 (ll/i -7^112-^7) , 
whence 

(/i,^.)2>^^^>0 (3.8) 

for rj sufficiently small. 

We now fix ?7 G (0, 1) small enough such that fl3^ holds true, let n > 1 and take ( G H'^{Q)nC'^{Tl) 
with C = 0: we aim at constructing a perturbation of v in Ai^is- To this end, we define 

7.(0 := -7^^ and ^,{t) := E [v + X{^,{0 + t) + X Q - [3 
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for A > and t E [—1/n, 1/n]. Since 

||(7.(C)+t)^. + C||2< + 11^,112 + 



(/i, z/^)2 n 
l(/i,C) 



< (1 + 2^^^) (l + ||^-7l||2) + ||Clh 



for all 77, > 1 and t G we have 

Mt) = - /3 + A (/i, IMO +t)u^ + C])2 + A £,(A) , 

where is a function which depends neither on n > 1 nor on t G [—1/n, 1/n] and satisfies e»)(A) — )■ 
as A — > 0. Owing to the property v G A^^^/j and the definition of 7,y(C)5 "we find 

( ±- ) = ±- (y", I'v)^ + A £r?(A) , 
\ n J n 

and it follows from the positivity (13. 8p of (/i, 2/^)2 that there exists A„ G (0, 1/n) such that 

V'A^ ( -- ) < < ( +- 
\ n J \ n 

The mean- value theorem then guarantees that there is t„ G (— 1/n, l/n) such that '^\„{tn) = 0. 
Setting C„ := (7r;(C) + ^n) i^r? + C) have thus shown that F[v + A„Cn] = /3- Since 



+ A„Cn = + A„(7^(C) + t„) I^r, + A„ C = « 



by (13.71) and the properties of (, we conclude that v + XnCn ^ ■^a,i3- Recalling that v, z/^, and ( belong 
to Hjf{Q), so that v + XnCn £ -^a,/3) ^^le minimizing property of u ensures that J^T,f[v] < ~^T,/[^+AnCn], 
that is 

\tzM + , E[.] < ^^-^'f-n'i ^rEl. + A„C.J . 
This inequality also reads 

< - {Av, W\v + A„,Cn) - W'iv))^ - A„ (AC, W'iv + A„,Cn) - ^v)^ 
+ ^ {W'iv + A„Cn) - W'iv), W'iv + A„Cn) + W'iv))^ 

+ ^ ||AC„||^ + ^ (i;-/,Cn)2 + ^ llCnIll (3.9) 
At It 

Observing that (Cn)n>i converges towards Coo := 7r?(C) i^r? + C i^i H'^i^) as n — )■ 00, we find that 

lim (ACn, W'iv + A„Cn) - Av) = [ fiAC^dx, 
lim A„ ||Cn||2 = lim A„ ||AC„||2 = 0, 

n— >oo n— >oo 
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while classical arguments and the embedding of H'^{Q) in L°°{Q) ensure that 

lim ^ {Av, W'{v + A„,Cn) - W'iv))^ = I W"{v) Coo At; dx , 

{W'{v + A„C„) - W'{v), W'{v + A„Cn) + W'{v))^ 

W'{v) W"{v) Coo dx . 



{V - /,Coo)^ 



r 
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lim — — 

n— >oo 2Xr 



Dividing fl3.9p by A,i and letting n ^ oo give 

< - /" W"{v) Coo^v dx - [ /i ACoo dx + I W\v) W"{v) Coo dx + 



fi [-A (7,(0 '^, + + W"{^) (7,(0 i^n + 0] dx + 



v-f 



T 



(7r?(C) + dx 



Observing that — C also belongs to Hjf{fl) fl C°°{Q) and satisfies (— C) = 0, the above inequality is 
also valid for — C and, since 7r;(— C) = ~'lv(0^ ^nd up with 



/i [-A (7,(0 i^r, + + W"{v) (7,(0 + 0] dx + 



(7,(0 iy, + dx = (3.10) 



for all C e H%{Q) n C°°{Q) satisfying C = 0. 

A first consequence of (13.101) is that /i G Hj^^Q): indeed, we can also write fl3.10p as 

/i AC c/x = - 7,(0 (/i, Ai/^)2 + 7,(0 (w"{v)fx + 
+ /iy»/i+^,C 



r 



Since v e H^{^), n e L^{^), and 



|7,(0I < 



l/^l|2 



it follows from the continuous embedding of H^{Q) in L°°{Q) that 



yU AC dx 



+ {\\W"(v 



oo ll/^lh 



W^"(i;) 
r 



1^112 



\v-fl 



< Civ,v,f,T) llClh, 

and a two-step duality argument entails first that fi G H'^{Q) and then that it satisfies the Neumann 
homogeneous boundary conditions, that is, 

^ieHl{9). (3.11) 
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We can then integrate twice by parts the first term of the left-hand side of (I3.10p to obtain 



[ E (7^(0 z/, + C) cia; = 
Jn 



for all C e Hl{n) n C^{n) satisfying C = 0, where 

S := - A/i + W"{v) /i . 

r 

We may now let 77 — )• with the help of (13. 7p and use a density argument to conclude that 

C - -112 (/^ - T^)) dx = 



for all C G L^(r2) satisfying C = 0- Alternatively, 

- A cdx=o. 

Now, take ( G L^{Q). Applying the above identity to C — C gives 



and, since this equality is valid for all test functions in L^{Q), we realize that 

^— ^ - A/i + W"{v) /i = S = A + 5/i in n, 

T 

with 

A:=E-^4^^^JI and 5 := f^^^ • (3.12) 

11/^-/^112 11/^-/^112 

The proof of Lemma [3.31 is complete. □ 

Noting that the Lagrange multipliers A and B arising in Lemma 13.31 are defined in a somewhat 
implicit way according to (I3.12p . the next step is to obtain some estimates on both of them. As we 
shall see now, this is quite easy for A + Bjt for which we have an explicit formula but turns out to 
be more complicated for B. 

Lemma 3.4. Assume that v G A^^^ solves the minimization problem (13. 2p and is such that fi : = 
— At> + W'{v) is not a constant. Assume further that f = a and Za,i3 = and consider a positive 
real number M > E{f). Then there is a positive real number km > depending on M such that 

\A + B-p\< Km and \B\ < km ( 1 + ^hjzlll\ . (3.I3) 
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Proof. Integrating ( 13. 5p over Q and using the properties fi G H^f^Q) and v = f = a, we obtain the 
identity 



A + BjL = W"{v) /i . 



It then follows from the Holder inequality, the continuous embedding of H^iVl) in L'*(n), f l2.ip . 
F{v) = /3, and (E3D that 



\A + B-[l\<C \\W"{v)h ||/i||2<C(l + ||t;||2) 
<C{l + \\v\\l,) ^/E[r]<CVM, 

which is the first inequality in f l3.13p . Next, since the average oi — fx over VL is zero, we can apply 
the operator M introduced in (12. 5p on fi — Ji and the function M {fi — Jl) belongs to Hj^^Q). We then 
infer from (12. 5p and (13. 5p that 

B\\VAf{fi-]l)\\l = -B I Af{fi--p) AAT (/i-/i) dx = B / (/i - /Z) A/" (// - /Z) dx 

Jn Ju 

= J (^^-^ - Afx + W"{v) fx- Ar(/i-7J) dx 

= j (^-^ + W'\v) 1^ A/'(/W-7i) dx+Wii-JlWl . (3.14) 

On the one hand, since Za^p = 0, f G A^^,^, and E[v] < E[f] < M by (13. 3p . we can apply Lemma [273] 
which simply yields 

\B\\\WAfil2-Jl)\\l>mM\B\. (3.15) 
On the other hand, according to (13. 3p . we have 

||/i-7l|l2 < 11/^112+ M\2 = 2||/w||2 = V^E[v] < C ^/E[f]< C \fM, (3.16) 
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which, together with Lemma [2.21 the Poincare-Wirtinger inequahty ( \2.3\} . the embedding of H^{Q) 
in L°°{n), ([23]), (Ea), and ([316]), gives 



v-f 



+ W"{v) 12] Af{l2-Jl) dx 



< 



< C2 

< C 

< c 

< c 

< C{M) 

We then infer from flXTI]) - flXT7]) that 

mM \B\ < \B\ ||VAr(^- 71)11^ 



+ c (1 + Ml) V^m] \\\/^^i^^-JI)\\ 





-/II2 


{ 


r 




-/II2 




T 




-fh 




T 




-fh 




T 




-fh 


\ T 



+ \\W"(v) 



+ (1 + 11^11^0 vm Wf^-J^ 



{1 + E[v]) VMj VM 
il + E[f]) ^^M] 



T 



(3.17) 



< 



v-f 



In \ 'T 
< C{M) 



+ W"{v) 12) Af{i2-Jl) dx 



which readily gives the second inequahty in f l3.13p since mjv/ > by Lemma 12.31 

4. Existence 



□ 



In this section, we prove the existence part of Theorem 11.11 We fix a G M and f3 G {f3a, 00) such 
that 

^a,/3 = 0. (4.1) 

Consider an initial condition t>o G A^^,/? ^^'^ ^ time step r G (0, 1). We define a sequence {v'^)n>i 
inductively as follows: 



^0 = ^0 , 

v^^i is a minimizer of Jv,i,j in A^^^ , n >0 , 



(4.2) 
(4.3) 



the functional J^t,vi being defined in (13. ip . Owing to Lemma 13. 1[ this sequence is well defined. 
Setting 



-Av: + w (v: 



n > 0, 



(4.4) 



A PHASE-FIELD APPROXIMATION OF THE WiLLMORE FLOW WITH VOLUME AND AREA CONSTRAINTS 15 

we define two piecewise constant time-dependent functions v'^ and /i"^ by 

K(t),/i-(t)):=«,/i;), tGK(n + l)r), n>0. (4.5) 
From tfie minimizing property (14.31) of f^+i, n > 0, we deduce the following estimates on {v'^,fi'^). 
Lemma 4.1. For r G (0, 1), ti > 0, and t2 > ti, we have 

Ey{t2)]<Ey{h)]<E[vo], (4.6) 

11^^(^1)11^.. + < Ci (l + , (4.7) 

oo 

Y,H+i-<\\l<2r E[vo], (4.8) 

n=0 

\\v^ih) - v^ih) \\l < 2E[vo] (r + t2 - h) . (4.9) 

Proof. Consider n > 0. According to tlie definition (14. 3 p of t'^+i, we have J^t,vi [^n+i] — •^t,vJ-^ [v^], 
that is, 

^%^ + i?K-M]<^K]- (4-10) 

On the one hand, the time monotonicity (14. 6 p oit E [t'^(t)] readily follows from (I4.10p by induction 
while the if^-estimate on v"^ and the L^-estimate on /i"^ in (14.70 are straightforward consequences of 
Lemma [2.21 the definition of E and /i'^, and (14. 6p . On the other hand, summing (I4.10p over n > 
gives fOj) . 

Finally, let ^2 > ^ and denote the largest integer smaller than tj/r by n^, z = 1, 2. We then 
infer from (14.80 that 



n2— 1 /n2— 1 



1/2 



n=ni 



whence (li^ . □ 

Consider next > 0. We observe that, since G H'^{Q), the function /i^+i defined in (14. 4p is 

constant if and only if —Av'^_^_^ + W — W' (t'^+i) = in f2, that is, v^^j^^ G -Z^^^. Since Z^^p 

is assumed to be empty, this situation cannot occur and we have thus established that /x^+i is not 
a constant for all n > 0. We are then in a position to apply Lemma 13.31 for each n > and deduce 
that jJin+i £ Hfjiyt) and there are real numbers Aj^+i and -B^+i such that /i^+i solves 



.J^±l-^-A^,l_^^ + W''{vl^,) ^^l^, = Al^, + Bl^,^il^, in fi. (4.11) 

We then define two piecewise constant time-dependent functions and B'^ by 

(A^(t),i?^(t)):=(A;,i?;^) , t G [nr,(n + l)r), n>l, (4.12) 
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and collect bounds for these functions in the next lemma. 



Lemma 4.2. For r G (0, 1) and T > t, we have 

[ [W{tt + \B-{t)f + dt < C,{T) . 



(4.13) 



Proof. Owing to (14. ip and the time monotonicity (14. 6 p of E [f^], the assumptions of Lemma [3.41 are 
satisfied with M = E[vo] for all n > and we obtain the estimates: 



\Al^, + B:^, < C and \B:^,\ < C \ 1 



n > 



(4.14) 



Let T > T and m > 1 be the largest integer smaller than T/r. On the one hand, we infer from (14.81) 
and fimj) that 



m— 1 



m— 1 



\B^{t)\' dt<T $^|5„V|'<Cr 5^ 1 + 



n=0 



n=0 



< C (mr + 2E[vo]) < C (1 + T) . 



(4.15) 



On the other hand, since v'^_^i G Hjf(il), we have /i^+i = W (f^+i) and, since v'^^i G it 
follows from the continuous embedding of H^{Q) in L^{Q) and Lemma [2.21 that 



< c 1 



^ n+l I 



Consequently, thanks to (I4.14p and (I4.15p . we have that 



m—l 



m—1 



n=0 



/ J IIL J. Ill, ± 

n=0 n=0 

<Cr + <C(1+T). 

Finally, we observe that 

/ ||A/i-(t)||^ rft<r $^||A^.^ 

^ r n 



(4.16) 



m—l 



''n+l\\2 ' 



(4.17) 



n=0 
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an , using ( I4.1ip and the continuous embedding of H'^{Q) in L°°(f2), that 













|2 
I2 












i2 


\K+i - 


n \ 


I2 












i2 






I2 



^ I II " '/i-t-J- "'^112 /-| II T 11^ \ ll ^ l|2 l/lT 1^ I D'^ |2 11 T 11^ 



hence, thanks to fl4.7p . we have 



i^i^UiW; <c\i + " + |A;^ir + i^n+ii 

for n G {0, . . . , m — 1}. We then infer from (14. Sp . f l4.15p . f l4.16p . fl4.17p . and the above estimate that 
£\\Afi^it)\\l dt<Cr (^m+'^^^p^^ + £^\\A^{t)f + \B^it)f) dt<C{T). 

Combining this estimate with fl4.15p . f l4.16p . and the L^-bound (14. 7p on /i^ gives (I4.13p . □ 

Thanks to the above analysis, all the tools required to perform the limit as r — are now available 
and we may thus proceed to identify the behaviour of (f as r — )■ 0. We begin with a consequence 
of (14 .Zp which guarantees compactness with respect to the space variable and (I4.9P which gives the 
time equicontinuity: owing to (14. 7p . (14. 9p . and the compactness of the embedding of H'^iVL) in H^iVL) 
and C (fi), a refined version of the Ascoli-Arzela theorem [H Proposition 3.3.1] ensures that there 
are a subsequence {v'^'')i^yi of {v"^) and a function 

veC ([0, oo) X n) n C([0, oo)- H\n)) n L~(0, oo; H\n)) (4.18) 

such that 

v^'it) — >v{t) in C{n)nH\n) forall t>0. (4.19) 

A straightforward consequence of (iTj), dun]), the continuity of the embedding of H^{n) in H^{n) 
and C (Q) , and Lebesgue's dominated convergence theorem is that 

t;"^ — yv in LP{0,T;C {Tl) Ci H\n)) forall pe[l,oo) and T>0. (4.20) 
In addition, it follows from (14. 7p and (I4.13P that we may assume that there are functions 

/i G L~(0, oo; L\n)) n L^O, oo; if2(f])) , ^ ^ ^^(0, oo) , i? G LL(0, oo) , 
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such that, for all T > (5 > 0, 

v^^^v in L°°{0,T-H\n)), (4.21) 

fi^" ^ 12 in L~(0,T;L2(fi)) and jj,^" ^ fi in L'^{5,T; H'^{n)) , (4.22) 

A^^ ^ A in L2(5,T) and B^^ ^ B in L'^{6,T) . (4.23) 

Now, since w^'^(t) G A^^^ and /i^'=(t) = -At;""''-(t) + VT' {v^^it)) for all t > 0, it readily follows from 
the convergences fHl^ . f lOT]) . and (KT2\i that 

t<t) G Ml^p and /i(t) = -At;(t) + W'{v{t)) for all t > . (4.24) 

It remains to derive the equation solved by f . To this end, we have to pass to the limit in fl4.1ip and in 
particular to identify the limits of the nonlinear terms {W" {v'^'') /i^'')fc>i and (i?'^*/i'^'=)^>^. Concerning 
the former, we combine the strong convergence (14.201) of iv'^'')i^yi with the weak convergence 04.221) 
of (/i'^'')fc>i to obtain that 

W" {v^') ^ W"{v) /i in L2((0, T) x fi) for all T > . (4.25) 

As for the latter, the situation is less clear as the convergences fl4.22p and fl4.23p are both weak 
convergences. However, we take advantage at this point of the fact that B'^'' depends only on 
time. Indeed, on the one hand, we notice that the strong convergence fl4.20p of {v'^'')k:>i implies 
that {fJ'^'')k>i converges strongly towards /i in Lp{0,T; H^^Q)') for all p G [l,oo) and T > 0. We 
combine this convergence with fl4.23p to obtain that [B'^'' fi'^'') converges towards Bfi in the sense 
of distributions. On the other hand, the sequence {B'^'^ jj^^) ^.^-^ is bounded in L^((5, T) x Vt) by (14. 7p 
and fl4.13p for all T > 5 > and is thus weakly compact in that space. Therefore, we have shown 
that, after possibly extracting a further subsequence, 

^r,. ^^fc ^ 5 ^ in ^2(^(^^ X fi) for all T > (5 > . (4.26) 

Now, for ^2 > ^1 > and G Iv^(f2), we denote the largest integer smaller than ti/Tk by rij^fc, z = 1, 2, 
and infer from fl4.1ip that 

n2,k 



•n.2,k „ 

n=ni,fe+l 
n2,k .(„+l)rfc I- 

/•(n2.fc+l)T"fc ^ 

= / / (A/i^'= - VT" (i;^'=) + A^^- + 5^^- /i^'=) (t) if dxdt 

J (rii t+Dri, J f2 
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Clearly, [rii^k + ti as A; — ?■ oo, i = 1,2. Letting /c — oo in the above identity gives, thanks to 

dEni), (|i22D, (1425]), g23D, and g26D, 



f {v{t2) - v{ti)) ipdx= r [ {Afi-W"{v) fi + A + B ij,){t) if dxdt. (4.27) 
Jq, Jti Jn 

It is now straightforward to check that Lemma I4l2] and the convergences (14.221) and (I4.23P imply that 

/i e L'^{0,T;H'^{n)), A G L'^{0,T), and 5 G ^^(CT) for all T > 0. Combining these integrability 

properties with (14.191) ensures that (I4.27P is also valid for ti = 0. The regularity of v and fi then 

allows us to deduce (fLTSj) from (14271) . 

It remains to check that A and B are given by (11.191) and (ll.20p . respectively: first, (I1.19P readily 

follows by integrating (ll.lSp and using that v(t) = a for all t > and the homogeneous Neumann 

boundary conditions satisfied by fi. Next, since -F[i;(t)] = /3 for all t > 0, we differentiate this identity 

with respect to time and, using once more the homogeneous Neumann boundary conditions for /i 

and the fact that 



fi {fi — fi) dx = 



n 



we obtain 



= / fidtv dx = - I W"{v) /i^ dx + {A + B-p) [ fi dx + B - Jl\\l . 

Jn Jn Jn 

The identity (ll.20p now follows from the above identity with the help of (I1.19p . 

Finally, fix t > and assume for contradiction that there is a sequence (s„)„>i in [0,t] such that 

II (/i — 7l)(s„)||2 — 7- as n — oo. Since [0,t] is compact, we may assume that s„ — > Soo as n — > oo 

for some Soo G [0,t]. Thanks to the regularity of v, we actually have v G C{[0,t]; H^{Q) fl C{Q)) 

so that ((yU — 'Jl){sn))n>i couvcrgcs towards (/i — /i)(soo) in H^{Qy. Since it also converges to zero 

in L'^{fl), we have shown that (/i — /i)(soo) = which implies that v{soo) G Za,i3 and contradicts 

(I1.17p . Therefore, — Ji\\2 is bounded from below by a positive constant in [0,t] and the proof of 

the existence part of Theorem 11.11 is complete. 

5. Uniqueness 

Consider a G M and (3 G (/3a, oo) satisfying (I1.17p . Let Vi, i = 1,2, be two solutions to (ll.7p -( |TTTT]) 
with fii := —Avi + W'{vi) and associated Lagrange multipliers {Ai,Bi), i = 1,2. Owing to the 
regularity of Vi and /ij, i = 1,2, stated in Theorem 11.11 and the embedding of H'^{Q) in L°°(f2), the 
function defined by 

(j){t) := |5i(t)| + \B2{t)\ + ||/ii(t)b2 + h2(t)\\H2 , t>0, 

satisfies 

0GL2(O,t) and ||/ii(t)||oo + ||/i2(i)||oo < (^0(t) for all t>0. (5.1) 
Also, given T > 0, Theorem ll.il (in particular (ll.2ip ) and the embedding of H'^{Q) in L°°(r2) ensure 
that there is Kj- > 1 such that 

||wi(t)||oo + lk2(t)||oo + ||/t^i(t)||2 + ||/U2(t)||2<i^T for all t G (0, T) , (5.2) 
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> — >0 for all tG(0,T) 



(5.3) 
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and 

mill 

Now, Vi — V2 solves 

dtiyi - V2) - A(/ii - /i2) = -W'iyi) III + W'\v2) /i2 + - + 5i /ii - B2 /i2 (5.4) 

in (0, 00) X Vt with homogeneous Neumann boundary conditions for vi — f 2 and ^1—^2- As a first step 
of the uniqueness proof, we estimate some terms in the right-hand side of (15. 4p . A first consequence 
of dED is that 

\\W"{V{) /ii - W"{V2) /i2|l2 < ||W^"(^'l)|loo 11/^1 - /^2||2 + \\W"{v,) - W"{v2)\\, 

< \\W"h^^^Kr,K,) ll/ii - /i2||2 + C0 - i;2||2 



<C(T)(||/ii-/i2||2 + 0||t;i-t;2|l2) 
It then readily follows from (15. 5p that 



iy"(t;i) /ii - W"{V2) /i2 < C(T) (ll/ii - ^211^ + - 

We next estimate Bi — B2. To this end, we first use (ll.20p to compute 

4 

\\\f^i-W\\l - 11^2 - A^lla B2\<^Ij, 



where 



A := |l|V/ii| 



2 - l|V/i2||2| , h 



\W"{vi) III - W"{v2) i4\ dx, 



h:= W"{vi) lii-W"{v2) li2 

Integrating by parts we find 

h < 



W"{V2) fl2 ||/il-/i2|h 



V(//i + /U2) ■ V(/ii - /i2) dx = {jii - /i2) A(/ii + /i2) c^a; 

<||Atl + At2||//2 ll/ii - /i2||2 < ll/il - /i2||2 • 

We next deduce from (15. 2p and (15.50 that 

h <||/il||2 ||l^"(t^l) /il - W"{V2) /i2||2 + ||W^"(t;2) /i2|l2 11/^1 ' ^^2\\2 

<C{T) (ll/ii - /i2||2 + - 1;2||2) + l|W^"(^'2)IL Il-"2||2 ll/il - y"2||2 

<C(T)(||/ii-/i2||2 + 0||i;i- 1^2112) • 
The last two terms are easier to estimate and we use (15. 2p and (15.61) to obtain that 

h + h <m'/^ ||/il||2 \W"{V,) /il - W"{V2) /X2I + ||l^"(t'2)IL U2\\2 11/^1 - /i2|h 
<C(r)(||/Xi-/i2||2 + 0|kl-^;2|l2) • 



(5.5) 
(5.6) 

(5.7) 



(5.^ 



(5.9) 



(5.10) 
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Using (IE2D and (lEl, we deduce from flO)-f l5lOD that 

i„ \\\fii-jH\\l Bi-\\fi2-M\l l^al I 2 II _||2| 

\Bl - B2\ <J —2 + —2 111/^1 - /illls - 11/^2 - M\2\ 

||/il-/il||2 ||/ii-/ii||2 

4 

<K^ + (p ll/il -7IT + /i2 -7i2|l2 ll/Ltl - /UT- /Lt2 + 7I2II2 

i=i 

<C(T) (1 + 0) (11^1 - ^2||2 + IK'l - V2\\2) ■ (5.11) 

After this preparation, we multiply (15. 4 p by vi — V2 and integrate over Q to obtain 

(t;i - V2) {W"{vx) /ii - W"{y2) /i2) dx 

n 

+ {Ai - A2) / {vi - 1^2) dx 



+ - i?2) / /ii (wi - W2) dx + B2 / (/ii - /i2) (vi - V2) dx . 
Ju Jn 

Since 

- ifii- fi2) A(t;i - V2) dx = {fii - fi2) (/ii - W^'(t^i) - /i2 + ^"'(^^2)) c^a; 

>\\f^l- f^2\\l - ll/il - /i2||2 ||W-'(t;i) - W'{V2)\\2 

>^ 11/^1-/^211^-^ 11^1-^211^ 

>\ ||/ii-/i2||^-C(T) ||t;i-^2||^ 
by (15. 2p and Young's inequality, we infer from vl = V2 = a that 
\ ^i2\\l -C{T) \\v,-V2\\l 

< \\Vl - t;2||2 \\W"{V,) /ii - V^"(t;2) /i2||2 + |5i - 52I Wfiih \\V, - V2\\2 
+ \B2\ ll/il - /i2||2 ||^'1 - ^'2||2 • 
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Using Young's inequality and <\5.2h . (I5.5p . and ( 15. lip , we deduce 

^11^1 -^2||2 + 11/^1 -/^2||2 

<C(T) \\v^-V2\\l + C{T) (1 + 0) (\\^,-^,\\^ + \\v,-V2\\2) \\Vi-V2\\2 
+ (j) ll/ii - /i2||2 IK'l - V2\\2 



1 



<2 ||/il-/i2||2 + C(^) (1+0') 11^1 -^^2 



Consequently, we obtain 

|(ll(^i-^2)(t)||^ + ^ £||(/ii-/i2)(.)il^rf.) <C(T) (l + 0^(t)) ||(^;i-t;2) 
and, since (1 + 0^) G -^^^(0, T) by (15.11) . the Gronwall lemma gives 

\\ivi-V2m\\l + l J\\i^^l-^^2)is)\\lds<CiT)\\iv^-V2m\\l for all tG[0,T]. (5.12) 
The uniqueness statement of Theorem 1 1 . 1 1 1 hen follows. 

Remark 5.1. By (I5.12p we have actually established a property of continuous dependence of the 
solutions to (ll.7p - (ll.lip on the initial data. 
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